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EVALUATION OF AN IMPEDANCE MODEL FOR PERFORATES INCLUDING THE 

EFFECT OF BIAS FLOW 


Juan F. Betts', Jesse I. Foiled, Jeffrey J. Kelly* 

Virginia Polytechnic Institute and State University 
Virginia Consortium of Engineering and Science Universities 
303 Butler Farm Road, Suite 101 
Hampton, VA 23666 

and 

Russell H. Thomas 5 
NASA Langley Research Center 
Hampton, VA 23681 

A new bias flow impedance model is developed for perforated plates from basic principles using as little 
empiricisms as possible. A quality experimental database was used to determine the predictive validity of 
the model. Results show that the model performs better for higher (15%) rather than lower (5%) percent 
open area (POA) samples. Based on the least squares ratio of numerical vs. experimental results, model 
predictions were on average within 20% and 30% for the higher and lower (POA), respectively. It is 
hypothesized on the work of other investigators that at lower POAs the higher fluid velocities in the 
perforate’s orifices start forming unsteady vortices, which is not accounted for in our model. The 
numerical model, in general also underpredicts the experiments. It is theorized that the actual acoustic C D 
is lower than the measured raylometer Cd used in the model. Using a larger Cd makes the numerical 
model predict lower impedances. The frequency domain model derived in this paper shows very good 
agreement with another model derived using a time domain approach. 


Nomenclature 

c speed of sound in the normal incidence 

tube 

Cd discharge coefficient 

d perforate hole diameter 

f frequency 

J 0 Bessel function of zeroth order 

J-i Bessel function of first order 

J 2 Bessel function of second order 

k 2itf/c 

L length of cavity 

M Mach number in the cavity 

M h effective Mach number in the holes of the 
perforate plates 
p pressure in the fluid 


p a acoustic pressure in cavity 

p a h acoustic pressure in perforate’s orifice 

t perforate thickness 

T period 

v fluid particle velocity 

v a acoustic particle velocity in the normal 
incidence tube 

Vah acoustic particle velocity in perforate’s 
orifice 

v b bias flow velocity in the cavity 

Vbh bias flow velocity in the perforate' s orifice 
v ms rms acoustic particle velocity 

X normalized reactance 

(i absolute viscosity of fluid 

\i' effective viscosity 


* Graduate Research Assistant, Mechanical Engineering Department, Member AIAA. 

f Graduate Research Assistant, Department of Mechanical, Aerospace, and Nuclear Engineering 
University of Virginia, Member AIAA. 

* Research Associate Professor, Mechanical Engineering Department, Senior Member AIAA. 

5 Aerospace Engineer, Aeroacoustics Branch, Senior Member AIAA. 

Copyright © 2000 by the American Institute of Aeronautics and Astronautics, Inc. No copyright is asserted 
in the United States under Title 17, U.S. Code. The U.S. Government has a royalty-free license to 
exercise all rights under the copyright claimed herein for Governmental Purposes. All other rights are 
reserved by the copyright owner. 


1 

American Institute of Aeronautics and Astronautics 



6 normalized resistance 

p density of fluid (air) 

a percent open area 

t time 

co 2jtf 

£, normalized impedance 

Throughout this paper acoustic impedance is 
normalized with respect to pc. The convection £=- 
Pa/pcvg and e*'“ T is used unless otherwise 
indicated. 


Introduction 

Due to the current and projected concerns 
about community noise annoyance from 
commercial air traffic, a number of enhanced 
abatement designs are being considered. One 
concept for improving suppression efficiency is the 
bias flow liner. Bias flow has been previously 
investigated by Dean, who conducted an initial 
proof-of-concept study that showed several 
potential advantages 1 . In-situ control of liner 
impedance would allow several desirable 
possibilities including optimizing liner impedance 
to match different operating conditions or to more 
accurately match design conditions that could not 
be reached due to manufacturing tolerances. 

The purpose of this study is to develop a 
bias flow impedance model from basic principles 
with as little empiricisms as possible and validate 
this model using quality experiments produced for 
this purpose from the NASA Langley Normal 
Incidence Tube (NIT). The experimental database 
included both fibermetal and perforate samples 
that were tested with and without bias flow in a 
single-degree-of-freedom liner configuration. 2 

A validation error criterion was developed 
and used as a standard for determining the validity 
of the numerical model results. 


accounted for by replacing the acoustic particle 
velocity in the impedance models with 


v . +v b (1) 

This correction seems more reasonable than 
Dean’s correction, especially at low bias flow 
velocities. Nevertheless, this bias flow correction 
as well as Dean’s correction does not follow from 
basic principles. 

In another study, Premo derived a bias 
flow model using a time-domain approach. This 
approach yielded the following replacement for the 
acoustic particle velocity to include bias flow 5 

A/( 1 -15v raB ) 2 +(2vJ ! (2) 

Development of A New Bias Flow Model 
The linear part of the following derivation 
follows from the principles outlined by Crandall of 
the impedance of a single tube. 6 The momentum 
equation for a viscous fluid is 


+ (v • V )v 

OX 


= -Vp + pV 2 v 


(3) 


where the compressibility term 


jpV(V.v) 


is 


assumed small and therefore is omitted from Eq. 
(3). Note, that by making this approximation 
density disturbances are ignored. Let any fluid 
variable q be the sum of a mean flow q b and an 
acoustic component q a , that is 


q=q b +q. (4) 

Substituting Eq. (4) into Eq. (3) and expanding 
yields 


Previous Bias Flow Corrections 
Dean developed a bias flow model in his 
initial proof-of-concept study (Ref. 1) of bias flow 
effect on liner impedance. His model is based on 
principles outlined by Hersh and Rogers. 3 Dean 
essentially replaced the acoustic particle velocity 
with the bias flow velocity in the Hersh and 
Rogers’ model without further justification except 
to say that the bias flow velocity was much greater 
than the acoustic particle velocity. As will be 
shown later, replacing the acoustic velocity with 
the bias flow velocity is questionable. 

In the NASA Langley Zwikker-Kosten 
Transmission Line Code (ZKTL) 4 , bias flow is 


^- + (v„»V)v b +(v b .V)v. 

[+(v. *V)v b +(v. *V)v. 

= -Vp„ -Vp. +pV 2 v b +(j.V J v i 


(5) 


This equation contains bias flow only, coupled bias 
flow/acoustic, and acoustic only components. The 
bias flow only components should balance, leaving 
coupled and acoustic only components. 
Furthermore, assume that perforate holes can be 
modeled as cylindrical ducts, and that the velocity 
has components in the x and r direction along the 
length and radius of the perforate hole 
respectively. So far, the above description pertains 
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to the interior of the perforate hole. To indicate 
this, the variables v a and v b will be subscripted 
with Vah and v bh . Making these approximations to 
Eq. (5) produces 


r dv . 9v h 9v 

p n*~ + v “ IT + v,b IT + v,h 


(6) 


.d£, 

dx 


A 

dr 2 r 9r dx 2 


Rewriting Eq. (6) and grouping linear and 
nonlinear coupling terms yields 


d x 


I <2* 


-H 


9V 


dr 2 


ldv^ + 3 2 v N 


dv 

+p l v “17 +v 

= _dp± 

dx 


r dr 

3v„ 


dx 2 


dv ^ 
+ v — - 
* dx * dx 


(7) 


The right hand side of Eq. (7) can be written as 

( 8 ) 


r 

ii 

s 

& 

fap* 


\ 

3x 1 

3x 

V 

k- dx 

nonlinear j 


Therefore the linear and nonlinear terms of Eq. (8) 
balance their corresponding counterparts 
respectively, in Eq. (7). By utilizing the definition of 
normalized impedance, the linear term of Eq. (8) 
can be rewritten in terms of the linear non- 
dimensional impedance by noting that 




dx 


K dv . 

“SltaorPC- — 


dx 


(9) 


Substituting the above relationship into Eq. (7) and 
assuming 3x=t, p ah =p a , dp a =Pa, and 3v a =v a 
produces 


-^i. a ,P CV *-p( V N, V 1 h)-^j :L = P 1 ( 10 ) 


P V . 


^P cv ! p 

C D 0 (C d o) ! V “ V - 2(C d o)’ 

= p. v . 


( 12 ) 


Integrating both sides of Eq. (12) over a period 
yields 


JLfi-PSU,.!}. 

Ti C„a TJ 


(C D o) ; 


■v.vjdi 


T J „2(C„o) j tJ 


(13) 


where v a and p a is assumed to be a harmonic 
solution with respect to time cos (cox). In Eq. (13) 
v’ is replaced by |v 4 |v, 2 so that this term remains 
positive. Simplifying this equation results in 


^tipeif P C v ? + 
C D a ~ 


P 2 

v k vl 


(C D a) : 

-2.828v^ =^pcv 


(14) 


2 

mis 


2(C„a) 2 3n 
Dividing both sides by pcv^ produces 




C D o 2 c(C d o) : 


[2v.+1.2v m ] 


(15) 


Melling notes that the pressure drop across a 
“sharp edge” orifice has been studied in some 
detail . From these studies, a departure from the 
nonlinear term in Eq. (15) is suggested. Making 
the corrections indicated by Melling to Eq. (15) 
yields 




l-o ! 


C D o 2 c(C d o) : 


-[2v b +1.2vJ 


(16) 


The linear impedance (linear) term is determined 
by solving Eq. (7) utilizing only the linear terms. 
Writing that equation leaving only the linear term 
yields 


From continuity, the relationship between the 
incident and hole velocities for an incompressible 
fluid is given by 


V , = C D CfV »h 


(11) 


Inserting the relationships in Eq. (11) into Eq. (10) 
and multiplying by v a produces 


9v.„ ( 9 ; v„ t 1 3v„ t 9'Vj, 

9r r 9r dx 2 




9p* 


(17) 


dx 


Melling summarizes Crandal’s solution to Eq. (17) 
assuming that the term d^'v^/dx 2 — » 0, which 
effectively eliminates any acoustic waves traveling 
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along the length of the duct. The impedance 
solution to Eq.(17) is 


_ _ 1 i27tft 
L 2 . 


(18) 


ao—1 .0 a^— “1.4092 32—0.0 

a 3 =0.33818 a4=0.0 a s =0.06793 

a^-0.02287 a 7 =0.003015 ae=-0.01614 

Note that the PBF model implies a bias 
flow correction to existing impedance models, 
where the acoustic particle velocity is replaced by 


where 

T (19) 

) 

Sivian showed that the viscosity inside the 
perforate hole along a highly thermally conductive 
wall was different from the absolute viscosity 
outside the hole 8 . It was shown in this study that 
for air over a wide range of temperatures, the 
relationship between p. and p' is 


k.d 


1 = 1 - 


k.-J„ 


M 

2 


2 v b + 1-2V™ (23) 

A problem arises when large enough negative bias 
flow rates would lead to negative impedance 
values. A recent study has experimentally shown 
the effect of negative bias flow rates of two 
degree-of-freedom liners. 11 The results of this 
suggest that increasing negative bias flow rates do 
not necessarily lead to lower or negative 
impedances. Consequently, Eq. (23) needs to be 
restricted to positive values. A possible method to 
accomplish this is to square Eq. (23) and take its 
square root. Doing this produces 


p' = 2.179p (20) 


v a =V 144v L+4v^+4.8v nm v b (24) 


Noting that the analysis for the impedance was 
done within the perforate’s hole ks becomes k'. 
Tijdeman has studied extensively the propagation 
of sound waves in cylindrical tubes 9 . He presents 
results of acoustic propagation in ducts under less 
restrictive assumptions, than what has been 
presented here. His results could be used as a 
starting point in developing a more comprehensive 
impedance model. 

Melling also describes contributions by 
Sivian and Ingard for end effects to the perforate 
hole, and Fok’s contribution to interaction effects 
between holes. This analysis is not going to be 
repeated here, but the results will be included. 
Therefore the Perforate Bias Flow (PBF) model is 




coC„ 


i27tft 


_ 2 _ 


3t cP 


8d 
k cT 


\V(a) 


( 21 ) 


+ (l r q M [ 2 v t + l,2vj 
2c(oC d ) 1 6 


and the Fok function \\i'{a) is defined as 10 


V (a) = ^a„(Vo) n (22) 

n=0 

where 


Thus the sign of v b does affect the magnitude of 
impedance but cannot drive it negative. Therefore, 
blowing produces a different impedance result 
than suction. 

Eq. (21) has a low and high frequency 
approximation. These approximations are valid if 


<i 27tfp 

2\ 


<1 


d 2rtfp 
2 \ 


>10 


(25) 


respectively. The low and high frequency 
approximations of Eq. (21) are 


t 32(it (l-c ! ) r 

t = — 4 — frl2v +1.2v 

PcoC D d 2c(oC D ) 


+ 1 - 


2nf 


acC n 


4 8d 

— t +- 


3 3th 


(26) 


£, = 2.82- 


r 



J 2d 

\y 

1 f 

V 1 

/ p 

1 t 


coC t d ^(oCj 


[2v h +\.2\ mm ] 


iIEL, 2 . 82 .fj /p) f l 

coC n caC. d 


8 d 


coC, d v'(a) 


(27) 


respectively. The PBF model runs into the problem 
that the Bessel functions in the linear component 
of the impedance obscure its physical relationship 
to geometrical parameters such as a, t, and d. The 
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low and high frequency approximations to the PBF 
model do not account for intermediate 
frequencies. Therefore, a single model that does 
not contain the Bessel functions, and “works” for 
all frequencies is desirable. One such correction 
similar to one done by Kraft, Yu, and Kwan, is 12 


Using the continuity principle for an 
incompressible fluid to relate the incident velocity 
and the hole velocity and Eqs. (29) and (30) yields 


Kl= 


I 

PJ0 SI% V2 

2coC d 

PC 


+ 2v b 


(31) 


16 jit 


- + 2.82- 


A j27d 

f f/1 


i/pj 


pcaC D d 2 caC d 

(i-+) 


2c(oC 0 )' 


-[2v b +1.2vJ + 


2icft 
ccrC ., 


+ 2.82 



12m */ 


2_i+ 8 d 


caC d d 3n y'(a) 


(28) 

Equation (28) will be called the Perforate Bias 
Flow Intermediate Frequency model or PBFIF 
model. 

Effective Mach Number 
M h is the effective Mach number at the 
perforate hole. It is defined as 


where |M H |is used to indicate the magnitude of 
the Mach number in the perforate. 

Validation Error Criteria 

The method of least squares can be 
employed to establish the validation criteria for the 
model predictions. Let E, and Nj be a set of 
experimental and numerical results, respectively. 
Assume that the following relationship exists 
between Ej and N| 13 

N, =SE, +e, (32) 

where ei is the random error associated with 
experiment i and S is an arbitrary constant. 
Minimizing the square of the errors yields the 
linear regression of N es timate on E where 


V + 2v 

M„= * — (29) 

The problem with using this equation is that the 
actual v^ is unknown as far as the impedance 
models are concerned. In fact that is what the 
impedance models are trying to find. The 
experimental Vah could be acquired from the 
experimental results and an assumed relationship 
between the incident and perforate hole acoustic 
particle velocities. Although this method would be 
the most accurate, a simplified approach suffices 
for this study. 

It is of interest to find out at what level of 
bias flow, do changes in M H due to SPL changes, 
not affect the overall acoustic impedance. 
Consequently, for a plane wave, the relationship 
between the acoustic particle velocity and 
pressure can be used to determine an 
approximate magnitude of v ah . Therefore, the 
relationship between the acoustic particle velocity 
and the SPL is 


P ,10 S% V2 


pc 


(30) 


S = 



(33) 


thus N es timate=SE j Note that this regression curve 
has been forced thru zero. S measures the ratio of 
numerical results to experimental results. If the 
ratio is equal to one then the experimental values 
equal the numerical results. Therefore, the 
deviation of this ratio from one is indicative of the 
least squares bias of the numerical vs. 
experimental results. 

The correlation coefficient is indicative of 
the random error in the predictions and is given 
by 14 


£(e,-eK-n) 

tel 

JS(E,-E)'i(N. -N? 

Y *=i 


(34) 


where, as before, E and N stand for experimental 
and numerical results, respectively. The 
correlation coefficient goes form -1 to +1 if the 
data is negatively or positively correlated, 
respectively. Note that in this case the regression 
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curve has not been constrained to pass thru 
N=E=0. 

Model Validation 

Figure 1 shows the experimental vs. the 
PBFIF model impedance results for a 
representative sample with no flow. The numerical 
results approach the experimental results as the 
sound pressure level is increased in both the 
resistance and reactance. The slope and 
correlation error criteria for both the resistance and 
reactance show continuous improvements as the 
sound pressure level is increased (See Table 1). 
This indicates that the nonlinear term in the PBFIF 
model, which is associated with the acoustic 
particle velocity better models the experiments, 
than the linear component of this model. 

Figure 2 shows the experimental and 
PBFIF model impedance at 130 dB and 5 % 
nominal (5.92 % actual) POA for various bias flow 
velocities. The figure shows increasing resistance 
with increasing bias flow rates. The resistance 
seems to stay constant with frequency except for 
the high bias flow rates, where it tends to 
decrease with increasing frequency. The 
reactance increases with frequency due to the 
cavity reactance. The experimental reactance 
tends to decrease as bias flow is increased. The 
model was not able to account for this effect. 

Figure 3 shows the slope and correlation 
error criteria for the resistance for various t/d 
samples. The flow rates are given in cm/s and M H 
(Hole Mach Number) for the NIT and hole 
velocities, respectively on the x-axis. The slope 
error criteria shows a mean slope starting at one 
for no flow then decreasing to 0.7 and increasing 
back to 0.9 as flow rate is increased. Therefore, 
the resistance numerical and experimental results 
are staying within 30% of each other in a least 
squares sense. 


SP-BxMUWIimn, M J906n, POA=5.9m> 


Experimental (Exp.) vs. PBFIF Model Numerical (Nun.) Results at 0 cm/s 




Frequency (ft) 


Figure 1. Experimental impedance vs. PBFIF 
model numerical predictions for a sample with 
POA=5.92. 


Table 1 . Slope and Correlation Coefficients results 
in Fig. 1. 





Resistance 



120 dB 

0.8127 

0.2509 

130 dB 

0.8692 

0.6413 

140 dB 

0.8670 

6.9290 

Reactance 



120 dB 

0.7552 

= 1.000 

130 dB 

0.7604 

= 1.000 

140 dB 

0.7865 

= 1.000 


The model underpredicts the resistances 
for all bias flow rates. At the higher flow rates the 
slope error improves because the experimental 
resistance decreases with increasing frequency at 
these flow rates (See Fig. 2). Since the model 
underpredicts the resistance, the decrease with 
frequency places the experimental resistance 
closer to the model results, and therefore 
decreases the measured slope error. 
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SP-52<d=0.65811 mm t t=0^90Gmm t tfd =1 .5052, POA=5.919%> 
Experimental (Exp.) vs. PBFIF Model Numerical (NumJ Results at 130 dB 

10 r 1 : ' T i 1 1 T 


0 
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« 
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3 
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1000 1200 1400 1600 1800 2000 2200 2400 2600 2800 3000 

Frequency (Hz) 

Figure 2. Experimental impedance vs. PBFIF 
model numerical predictions for a sample with 
POA=5.92. 

The mean correlation increases from 0.45 
to 0.65 and then goes negative. Part of the reason 
for the low correlation numbers is the low 
variability with frequency that these impedance 
curves have (See Fig. 2). Since the correlation 
measures the percentage change in experiments 
accounted by the model, low variation increases 
the errors. 

Figure 4 shows the slope and correlation 
error criteria for the reactance of the same t/d 
samples. The mean reactance slope starts at 0.7 
goes to 1 .1 and then decreases to negative values 
as bias flow is increased. The mean correlation on 
the other hand stays close to one at the lower bias 
flow rates and then sharply decreases to 0.65. 

Note that there is significant scatter in Fig. 
3 and 4 for the various t/d samples. This scatter is 
low at no flow, but significant for nearly all the flow 
rates. This result indicates that there is some sort 
of interaction between t/d and bias flow velocities 
not accounted for by the PBFIF model. 

Jing and Sun 15,16 provide possible 
explanations for the scatter in error associated 
with various t/d, and the impedance values at 
higher bias flow rates. They developed two 
models, one that looked into the nonlinear 


- Bp.OcnVs 
| -6- Eip. 100 cm/s 

* Bp. 200 cm/s 
j-$- Bp.30GcnV$ 

I -9- Bp. 600 cm/s 

• Mum. 0 cnVs 
| -o* Num. 100 cm/s 

- Niri200on/s 
I Nun.300cnVs 

-3- Nmv600cm/s 



u u - o-9- 


-- &—B— e--n e--a- - -o— c* ~g~s~ ■&— e— e— b— e— e— a— b~b— t> -i 


♦ - t - ♦ > 4 ♦ ♦ 4 ♦ 4 - 4 -4*4 M M 1" 


0 0 0 




properties of the hole and another that accounted 
for bias flow. 




Figure 3. PBFIF model resistance error criteria for 
POA=5. 


Jing and Sun (Ref. 18) took issue with the 
empirical nature of the determination of C D . Their 
model showed that the formation of the vena 
contracta is a time dependent problem in the 
presence of an acoustic field, and determined an 
approximate steady state value of 0.61 for Cd- 
This value of C D is lower than the normally 
measured experimental value of 0.76 or higher. 
This possibly explains why the PBFIF model 
underpredicts the resistance experiments since a 
higher C D produces lower predicted resistance 
results. 

Jing and Sun (Ref. 19) developed a 
boundary element bias flow model from Euler’s 
Equation, where the Kutta condition was applied to 
generate unsteady vortices due to the bias flow. 
Their model showed that different t/d and hole 
Mach numbers had a significant impact in the 
resistance and reactance. Their model 
successfully predicted a decrease in reactance as 
bias flow was increased. Since the PBFIF model 
does not take into account vortex creation in the 
perforate’s holes, it fails to predict the decrease in 
reactance, as well as, some of the interaction 
between hole Mach number and changes in t/d. 
This possibly explains the scatter seen in the error 
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plots with t/d, which suggest that the models are 
not accounting for some physical phenomena 
associated with t/d. 


Reactance Error Crtoia vs. Bes Fkiw Rate for P0A=5% 
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correlation between the experimental and 
numerical results. 

SP-54<d=0.63906inin, HMm, W =1.5898, POA»16.481%> 

Experimental (Exp.) vs. PBRF Model Numericai (Mum.) ResuBs at 1 30 dB 
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Figure 4. PBFIF model reactance error criteria for 
POA=5. 

The PBFIF model improves its prediction 
as POA is increased. Figure 5 shows the 
experimental and PBFIF model impedance at 130 
dB and 15% nominal (16.48% actual) POA for 
various bias flow velocities. The figure shows the 
experimental and numerical results are closer at 
15% than at 5 % (Fig. 2) POA. The mean slope 
and correlation error criteria for both the resistance 
and reactance show significant improvements as 
seen in Figs. 6 and 7. The error plots also show 
significant reduction in scatter due to different t/d, 
indicating that the PBFIF model can better account 
for changes in t/d at this higher POA. 

For the resistance, the mean slope criteria 
goes from about 0.9 at zero flow to 1 .2 and down 
to .8 in Fig. 6, indicating that the numerical and 
experimental results stay on average within 20% 
of each other. The correlation shows 

improvement, but still exhibits poor correlation in 
part due to the low variability in resistance curves. 
The reactance slope criteria starts around 0.8 and 
steadily increases to one, indicating that the 
reactance curves are within 20% of the 
experiments. The reactance correlation is about 
one for all flow rates indicating a nearly perfect 


Figure 5. Experimental impedance vs. PBFIF 
model numerical predictions for a sample with 
POA=16.48. 

There is a possible explanation for the 
PBFIF model’s improved performance as POA is 
increased. As POA is increased, the effective bias 
flow Mach number in the perforate hole 
decreases. For example, the highest Mach 
number achieved with the 15% and 5% POA 
samples was 0.126 (See Figs. 6 and 7) and 0.357 
(See Figs. 3 and 4), respectively. The PBFIF 
model performs better at lower flow velocities, 
because the formation of unsteady vortices is a 
function of the fluid velocity, given constant 
viscosity and hole diameter. At low flow velocities 
vortices are not formed, and the perforate’s hole 
physics follow the principles outlined in the PBFIF 
model. As the bias flow velocity in the hole 
increases, unsteady vortices start forming and the 
assumptions and principles outlined in the PBFIF 
model start to break down. At higher perforate bias 
flow rates a model such as that outlined in Ref. 19 
should be used. 
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Resistance Error Criteria vs. Bias Flow Rate for POA=15\ Reactance Error Criteria vs. Bias Flow Rate for POA=1 5X 




Bias Rcw Rate (cnVs^ 



Bias Flo* Rale (anteJKM^} 


Figure 6. PBFIF model resistance error criteria for 
POA=15. 

The PBFIF model implies a bias flow 
correction given by Eq. (23) as indicated earlier. 
This correction, which was derived using a 
frequency domain approach, is similar to that 
developed by Premo (Ref. 6) using a time domain 
approach (See Eq. (2)). Figures 8 and 9 show 
resistance comparisons between the models for a 
5% and 15% POA sample perforates. The 
reactance is not shown because bias flow only 
affects the resistance in both Premo’s and the 
PBFIF model. These figures show almost identical 
resistance values for all bias flow rates. 


Figure 7. PBFIF model reactance error criteria for 
POA=15. 

$P-52<d=0.S5811mm, t=0^906men, tld • =1.5052, POA=5.919%> 

PBFIF vs. Prow's Bias Flew Correction at 130 dB 
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Figure 8. PBFIF vs. Premo’s impedance model 
predictions at POA=5.92. 
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PBFIF vs. PremeTs Bias Flow Correction at 1 30 dB 



Frequency (Hz) 


Figure 9. PBFIF vs. Premo's impedance model 
predictions at POA=16.48. 

Conclusions 

The numerical results in general 
underpredict the experiments. The validation error 
criteria determined that the model performed 
better for higher (15%) rather than lower (5%) 
percent open area (POA) samples. The slope 
error criteria indicated that the model’s predictions 
were on average within 20% and 30% for the 
higher and lower (POA), respectively. The lower 
POA perforate numerical model predictions 
presented more scatter for different t/ds, in the 
slope and correlation error criteria, than at higher 
POA. The scatter also tended to increase as the 
flow velocity was increased for a given POA. 

It is hypothesized that the model 
underpredicts the experiments because the 
measured raylometer C D is higher than the actual 
acoustic C D as indicated in Ref. 18. Also, at lower 
POAs the higher fluid velocities in the perforate’s 
orifices start forming unsteady vortices, which is 
not accounted for in our model. Consequently the 
model performs better at the higher rather than 
lower POAs. The reader is referred to Ref. 1 9 if a 
higher bias flow velocity model is desired. Note 
Ref. 19 utilizes the e’ 10 ” while this paper utilized the 
e +kBt convention. 


The PBFIF model showed very good 
agreement with Premo’s time domain model (Ref. 
6) for all bias flow rates. 
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